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We report on an on-going two-flavor full QCD study on CP-PACS using an RG-improved gauge action and a 
tadpole-improved SW quark action. Runs are made for three lattice spacings a~ ps 0.9, 1.3, and 2.5 GeV on 
12 3 x 24, 16 3 x 32, and 24 3 x 48 lattices. Four sea quark masses having mps/mv ~ 0.8-0.6 are simulated, for each 
of which hadron masses are evaluated for valence quark masses corresponding to mps/mv « 0.8-0.5. Results for 
hadron and light quark masses are presented and compared with those obtained in quenched QCD. 



1. Introduction 

Having observed an unambiguous deviation of 
the quenched light hadron spectrum from the ex- 
periment |l|-|j| , we have started, as a logical next 
step, a systematic full QCD simulation. Due to 
the necessity of working at coarse lattice spacings 
to cope with an increased demand on computer 
power, we employ an RG-improved gauge action 
combined with a meanheld-improved SW quark 
action. This choice is an outcome of a compara- 
tive study of various action combinations carried 
out prior to the present workQ. In this report 
we describe the main points of results obtained 
so far, referring to Ref. || for details. 

2. Parameters of simulation 

We study QCD with two flavors of sea quarks, 
identified with the degenerate u and d quarks, 
treating the strange quark in the quenched ap- 
proximation. Based on test simulations, we 
choose the run parameters listed in Table [l] cor- 
responding to three lattice spacings in the range 
a- 1 ~ 1-2 GeV. We employ 12 3 x 24, 16 3 x 32, 
and 24 3 x 48 lattices to keep the spatial size at 

*prcsented by K. Kanaya at Lattice98, Boulder, Colorado, 
USA, 13-18 July 1998. 



Table 1 

Parameters of simulation, 

lattice a -1 [GeV] K sea m PS /m v N con{ 

P/c SW La[fm] XjVsepr 

12 3 x 24 0.917(10) 0.1409 0.806(1) 1250x5 

1.80/1.60 2.58(3) 0.1430 0.753(1) 1000x5 

0.1445 0.696(2) 1400x5 

0.1464 0.548(4) 1050x5 

16 3 x32 1.288(15) 0.1375 0.805(1) 1400x5 

1.95/1.53 2.45(3) 0.1390 0.751(1) 1400x5 

0.1400 0.688(1) 1400x5 

0.1410 0.586(3) 1000x5 

24 3 x 48 2.45(9) 0.1351 0.800(2) 250x5 
2.20/1.44 1.93(7) 0.1358 0.752(3) 270x5 
0.1363 0.702(3) 322x5 
0.1368 0.637(6) 253x5 



L ~ 2.4 frn. Simulations are made at sea quark 
masses corresponding to mps/my « 0.8, 0.75, 
0.7 and 0.6. 

At each sea quark mass, the light hadron spec- 
trum is computed using valence quark masses cor- 
responding to mps/my ~ 0.8, 0.75, 0.7, 0.6 and 
0.5. We smear the quark source with an expo- 
nential smearing function Q. All unequal quark 
mass combinations allowed for degenerate u and 
d sea quarks are taken. Hadron masses are ex- 
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Figure 1. Vector meson mass at j3 — 1.95. Lines 
are results of quadratic fit. S and V correspond 
to K va i = K sea and ^ K sea , respectively. 

tracted by uncorrelated fits. Errors are estimated 
by the jackknife method with a bin size of 10 con- 
figurations generated over 50 HMC trajectories. 

3. Hadron spectrum 

In Fig. Q we plot the vector meson mass at 
/3 = 1.95 as a function of the average of 1/K va i 
of the valence quark pair of the meson. Open 
points along each dashed lines, which are re- 
sults of global fits described below, are "partially 
quenched" results calculated for five values of 
K va i for a fixed K sea of sea quark. While each of 
these mass results are almost linear, their slope 
show a variation with K sea . Moreover, the full 
QCD results satisfying K va i = K sea , which are 
the four filled points, exhibit a concave curva- 
ture. Similar behavior is observed also in other 
hadrons. 

We then make a fit of the mass results with a 
general quadratic ansatz in 1/K sea and l/K va [. 
For pseudoscalar mesons the fit ansatz is 

mp S = B s fh sea + B v fn va i + C s rh 2 sea + C v m 2 val 

+C sv m sea m va i + Ci2rh va i(iyn val ( 2 ) 

where bare quark masses are defined as 
m sea/vam = {K^ a/val{l) - K~ 1 )/2 with m va i = 
(m„ a ;(!) +m va i( 2 ))/2 the average valence quark 
mass. Similar Ansatze without the C12 term are 
used for vector mesons and decuplet baryons. For 
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Figure 2. Hadron masses as a function of lattice 
spacing. Open squares are the results of quenched 
QCD with Wilson action @. 

octet baryons a form inspired from chiral pertur- 
bation theory with general quadratic terms of in- 
dividual m va i(i) are employed, simultaneously fit- 
ting A and E-like baryons. 

From the fit results, we determine the physi- 
cal light quark point K u d from m w /m p , and the 
strange quark point K strange from rriK, ra^*, or 
m^. The scale is set by m p at K u( i. Physical 
hadron masses as a function of the scale are shown 
in Fig. ^ where the results of quenched QCD with 
the Wilson quark action are also plotted (open 
symbols). 

A very interesting indication in the meson sec- 
tor is that the present two-flavor full QCD re- 
sult for the hyperfine splitting extrapolates to 
a value noticeably closer to experiment than 
that for quenched QCD in the continuum limit. 
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Figure 3. Light quark mass computed with vari- 
ous definitions. 



The remaining discrepancy might be due to the 
quenched treatment of the strange quark itself. 

Implications from the baryon results are less 
clear. While scaling violations are small, the mass 
results lie 5-10% high compared to experiment. 
We note two points in this regard : (i) the lattice 
size of 2.4 fm (actually turned out to be 1.93 fm 
at f3 = 2.2) is smaller compared to 3.0 fm for 
the quenched case, hence finite-size effects can 
be an issue, and (ii) the smallest quark mass for 
full QCD corresponds to mps/my ~ 0.5 rather 
than w 0.4 for quenched QCD, possibly leading 
to an overestimate of baryon masses. Further 
work matching these points is needed to clarify 
sea quark effects in the baryon mass spectrum. 

4. Quark masses 

Quark mass can be defined either by vector 



Ward identity (VWI) m\ 



VWI 



i , with m q = 

(1/K q — l/K crit )/2, or by axial Ward identity 
(AWI) ZaV^ = 2mf Nl Z P P. A recent discus- 
sion with the VWI definition has been whether 
to take the critical value K s c ea of sea quark for 
K cr it or the partially quenched critical value 
defined by m PS (K va i) = with K sea = K ud 
fixed The former is a natural choice for the 
averaged mass m U d of u and d quarks since sea 
quarks in two-flavor full QCD are identified with 
them. The SESAM Collaboration found, how- 
ever, that this choice leads to a large value of the 
ratio m s /m u( i ~ 50 at a -1 w 2 GeV. The value 



of m u d itself differs by a factor two depending on 
the two choices. This effect is negiligibly small 
for the heavier strange quark. 

In Fig. H we show our results for m v d at 2 GeV 
in the MS scheme. Renormalization coefficients 
are taken from one-loop tadpole- improved pertur- 
bation theory ||. The results for the two choices 
of K cr n in VWI, while sizably differing at finite 
lattice spacings, converge to a common value of 
2.5-3 MeV in the continuum limit, and so does 
the AWI quark mass. This value is about 40% 
smaller than 4.6(2) MeV from quenched QCD [§. 

The convergence found above also resolves the 
problem with the ratio m s /m u d for the choice 
K cr it = Kc ea ; while we find m s /m u d ~ 50 at 
a -1 ~ 0.8 GeV, the ratio decreases with the lat- 
tice spacing, converging in the continuum limit 
towards a value m s /m u d ~ 25. 

For the strange quark mass itself, the main un- 
certainty arises from the choice of input to fix 
K strange- We find m s w 70 (m^-input) - 80 
(m^ -input) MeV in the continuum limit. Com- 
pared to 115(2) (rriK- input) - 143(6) (m^-input) 
obtained in quenched QCD(|], the discrepancy, 
which might be explained by the quenching of 
strange quark, is smaller, and the values are also 
smaller. 
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